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ABSTRACT 


For  a  communication  system  using  satellites  carrying  antennas  tied  to  a  switching  system, 
the  probability  distribution  of  angles  between  the  pointing  directions  of  the  axis  of  the 
operating  antenna  and  of  a  ground  receiving  station  is  found.  This  study  is  based  on  the 
condition  that  the  antenna  pointing  most  nearly  in  the  direction  of  a  sensing  signal  is 
switched  on.  In  the  determination  of  the  characteristics  of  this  system,  geometries  of 
pointing  directions  are  found,  which  are  optimum  or  near  optimum,  for  the  cases  of  4, 
6,  8,  12,  16,  24,  32,  48,  and  60  antennas.  Consideration  is  also  given  to  determining 
an  optimum  antenna  pattern  to  be  used  with  any  geometry.  For  any  pattern,  a  conver¬ 
sion  from  probability  distribution  of  angles  to  probability  distribution  of  gains  is  straight¬ 
forward.  The  results  of  this  study  are  summarized  in  a  chart. 
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PROBABILITY  DISTRIBUTION  OF  ANTENNA  GAIN  FOR  SATELLITE 
WITH  SWITCHED  ANTENNA  SYSTEM 


I.  INTRODUCTION 

A  communication  system  is  presently  being  studied  that  employs  a  number  of  active 
satellites  dispersed  in  separate  near-synchronous  orbits  about  the  earth.  Each  satellite  will 
have  an  array  of  antennas  only  one  of  which  will  radiate  at  any  one  time  during  operation.  The 
satellite  must  carry  instrumentation  to  monitor  a  signal  from  the  earth  and  switch  on  the  antenna 
that  points  most  nearly  in  the  direction  of  the  signal.  On  the  other  hand,  the  need  of  satellite 
stabilization  is  eliminated.  Without  describing  equipment  design,  this  report  discusses  arrange¬ 
ments  of  various  numbers  of  antenna  pointing  directions  and  the  antenna  gains  in  the  direction  of 
an  earth  receiving  station.  The  RF  power  that  must  be  generated  in  the  satellite  should  be  as 
low  as  possible,  whereas  the  power  transmitted  from  a  ground  station  can  be  sufficiently  high  so 
that  the  location  of  the  transmitter  may  be  eliminated  as  a  factor. 

The  direction-sensing  system  that  provides  the  logic  for  determining  which  antenna  to  switch 
on  may  be  designed  to  monitor  the  signal  from  the  transmitter,  a  signal  sent  from  the  receiving 
station  or,  if  an  infrared  system  is  used,  the  signal  arriving  apparently  from  the  center  of  the 
earth.  This  direction  will  be  denoted  as  OS  with  the  satellite  situated  at  O,  as  shown  in  Fig.  1. 

The  satellite  carries  N  antennas  all  of  which  will  be  assumed  to  have  the  same  circular 
symmetric  antenna  pattern.  The  pointing  directions  of  the  centers  of  their  patterns  will  be  de¬ 
noted  by  OT^  (n  1 , 2, .  .  .  ,  N)  and  that  of  the  antenna  in  operation  by  OT. 


Fig.  1.  Schematic  view  of  communication  system  showing  operation  of  one 
satellite  O.  The  signal  sensed  propagates  along  direction  OS  and  the  an¬ 
tenna  whose  pointing  direction  OT  is  most  nearly  in  direction  OS  is  switched 
on.  The  receiver  is  in  the  direction  OR. 
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Other  designations  will  be 

OH  (direction  of  the  receiving  station), 
o?  (angle  between  OR  and  OT), 

0  (angle  between  OH  and  OS), 
y  (angle  between  OT  and  OS). 

In  normal  operation,  it  can  be  assumed  that  of  the  directions  OT^,  OT  subtends  the  smallest 
angle  with  OS,  that  is,  the  antenna  pointing  most  nearly  in  the  direction  of  the  signal  sensed  will 
be  the  one  in  operation.  Hence,  the  region  enclosing  the  satellite  can  be  pictured  as  divided  into 
pyramidal  segments  by  planes,  each  of  which  bisects  the  angle  between  two  adjacent  directions 
OT^  and  is  normal  to  the  plane  in  which  the  two  OT  lie.  Each  segment  defines  the  region  in 
which  OS  must  lie  for  the  corresponding  antenna  to  operate. 

It  is  apparent  that  for  any  geometry  of  the  pointing  directions  OT  ,  angle  y  has  a  maximum 
value  vm.  P'urthermore,  the  maximum  value  of  angle  a  would  be  p  +  vm-  In  mathematical 
terms 

P(a  <  a)  =  1  ,  a  >  /3  +  vm  (la) 

or 

P(cv  <  a)  <  1  ,  a  <  p  +  vm  (lb) 

One  would  want  to  choose  the  geometry  that  gives  the  lowest  values  of  v^.  This  allows  the  an¬ 
tennas  to  be  designed  for  narrower  beamwidths  and  correspondingly  higher  gains.  The  HF  power 

source  requirements  would  then  be  lessened. 

The  concept  of  the  system  and  its  operation  has  now  been  made  clear.  To  resolve  the  value 
of  the  system,  the  discussion  will  follow  the  procedure  outlined  below: 

(a)  Development  of  the  general  equations  to  determine  the  probability 
distribution  P(a?  <  a),  the  probability  that  angle  cx  is  less  than 
various  angles  a  for  various  values  of  angle  0 . 

(b)  Programming  the  equations  for  calculation  by  a  computer. 

(c)  Selection  of  pointing  direction  geometries  for  different  numbers 

of  antennas  that  will  give  minimum  or  near-minimum  values  of  v  . 

The  cases  that  are  explored  are  N  =  4,  6,  8,  12,  16,  24,  32,  48, m 
and  60  antennas. 

(d)  Calculation  and  plot  of  P(a  <  a)  for  these  geometries. 

(e)  Determination  of  the  optimum  antenna  for  these  geometries  based 
on  a  representative  antenna  pattern  g  =  2J^(g0  sina)/sina. 

(f)  Consideration  of  the  HF  circuits  supplying  the  antennas  and  estima¬ 
tion  of  the  circuit  losses. 

The  results  of  this  discussion  are  summarized  in  Figs.  6  to  24,  Fig.  27,  and  Appendix  A. 

II.  DEVELOPMENT  OF  THE  EQUATIONS  FOR  P(ct  <  a) 

Consider  an  arbitrary  pyramidal  segment  (Fig.  2).  Introduce  an  incremental  sheet  bound  by 
two  cones,  TOS°  =  y  and  TOS°  =  y  +  dy,  and  the  boundary  planes  of  the  segment.  The  probability 
that  OS  will  fall  within  the  strip  is 
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1 3~6?-?595 1 


Fig.  2.  Geometry  associated  with  antenna  in  operation. 
It  is  characterized  by  a  segment  of  a  sphere  whose  bound¬ 
ary  is  made  up  of  planes  bisecting  the  angles  between  ad¬ 
jacent  pointing  directions. 


where  S2  is  the  solid  angle  enclosed  and  4ir  is  the  solid  angle  of  the  whole  celestial  sphere.  (This 
is  based  on  the  condition  that  all  directions  OS  are  equally  probable.)  In  turn,  the  solid  angle  £2 
is  the  product  of  the  angular  length  £  and  angular  width  dy.  Hence, 


P 

s 


*tdy 

47T 


(3) 


Angle  ft  is  measured  on  the  surface  of  the  cone  TOS°  =  y  at  its  apex  (as  if  the  cone  surface  were 
slit  along  a  generator  line  and  unrolled  before  the  measurement).  This  surface  is  sliced  by  the 
boundary  of  the  pyramidal  segment,  and  £  is  the  sum  of  the  angles  of  the  surface  sectors  inside 
the  boundary. 

For  any  direction  OS,  direction  OR  will  fall  on  the  surface  of  another  cone  whose  axis  is 
OS,  since  the  angle  ROS°  =  0  is  invariant.  (Each  value  angle  0  assumes  may  be  viewed  as  an 
independent  situation.)  The  probability  that  a  <  a  is 

£ 

P  =  ,  (4) 

a  £ 


where  £  is  the  angle  of  the  sector  measured  on  the  surface  of  the  cone  ROS°  =  0  for  which  a  <  a 
and  ic  is  the  total  angle  of  the  cone  surface  measured  at  its  apex.  (This  is  based  on  the  condition 
that  all  directions  OR,  where  ROS°  =  0,  are  equally  probable.) 

The  probability  that  direction  OS  lies  within  the  incremental  sheet  and  a  <  a  is 

£.£  dy 

p'  =  p  p  =  4A-  • 

s  a  4t r£ 

c 


If  direction  OS  fell  within  any  other  incremental  sheet,  the  probability  would  be  given  by  the 
same  formula,  the  values  of  £  and  £g  being  different.  The  probability  that  a  <  a  for  any  direc¬ 
tion  OS  must  then  be  the  sum  of  the  individual  probabilities,  or 

£  £  dy 

P<«<a)=Z^7-  • 

c 

The  summation  over  one  pyramidal  segment  is,  in  the  limit,  an  integration.  The  sum  over  all 
segments  is  a  sum  of  the  integrals.  So 
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P(o>  <  a)  = 


(6) 


£  £  dy 

y  \  _Lg_ 

L  J  4*1 
n- 1 

where  n  denotes  the  segment.  The  summation  can  be  carried  out  over  all  the  segments  for 
the  same  angle  y  before  integrating.  Angle  £^  would  be  the  only  factor  that  could  change  from 
one  segment  to  another.  The  formula  can  then  be  given  finally  as 


P(o  <  a) 


(7) 


It  remains  to  find  expressions  for  £  /£^  and  £ 

For  £g/£c,  Fig.  3(a-f)  shows  the  various  possible  situations  that  may  exist  and  the  value  of 
£g/£c  which  corresponds  to  each  situation.  Four  of  the  cases  are  self  evident;  the  other  two, 
a-/3<y<a+/3,  a>/3  and  fi  —  a  <  y  <  a  +  (3 ,  a  (3  are  derived  in  Appendix  B. 


a  >  p. 


(0  a-/9  <  y<Q  +/9 

ic  |  cos  a  -cos  b  cos  y 

«7=^C0S  sin  /9  sin  y 


~p- 62-2596 1 

a  <$: 


(d)  /9-a<y<a+/9 

is  i  *  cos  a  -cos  s  cosy 

— £  =  -L  fOC  *  1  - — _ L. 

*c  w  sin/9  sin y 


Fig.  3.  Various  situations  that  may  exist  for  the  relative  values 
of  the  angles  a,  p,  and  y.  The  value  £s/£  is  given  in  each  case. 
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Fig.  4.  A  detailed  portion  of  the  geometry  of  Fig.  2  for  the  purpose 

of  determining  f  . 

tn 


To  derive  the  equations  for  a  geometry  for  one  segment  is  set  up  as  shown  in  Fig.  4. 

The  segment  belongs  to  direction  OT,  and  its  boundary  planes  are  determined  by  adjacent  direc¬ 
tions  OT.  (i  =  1,  2,  .  .  .  ,  I)  with  i  increasing  as  one  goes  counterclockwise  around  the  line  OT. 
i  th 

Thus,  in  Fig.  4,  j  =  i  +  1,  or  if  i  =  I,  j  =  1.  The  i  boundary  plane  bisects  the  angle  between  the 
lines  OT  and  OT.  and  intersects  the  plane  containing  the  lines  at  right  angle  along  the  line  OF.. 

Furthermore,  adjacent  boundary  planes  intersect  along  lines  OG...  The  cone  defined  by  TOS°  =  y 
t  h  ^ 

may  intersect  the  i  boundary  plane  along  lines  OH^  and  OfT^.  The  subscripts  a  and  b  will 

be  used  to  denote  the  one  or  the  other  side  of  the  line  OF..  Let  points  T,  T.,  T.,  F.,  F  ,  G.  , 

l  ^  i  J  1  J  ij 

H.^,  and  11.^  all  lie  on  a  plane  which  is  perpendicular  to  direction  OT.  In  addition,  there  are  the 

designations  TOF.°  =  u  .,  TOG..°  =  v  ..,  F.TH.  °  =  (p  .  ,  F.TG..°  =  0  ..  ,  as  shown  in  Fig.  4, 
l  ni'  ij  nij  l  la  rma  i  ij  mb' 

and  T.OT . °  =  2u  ...  Also,  TF.G..°  =  tt/2. 
i  J  nij  i  ij 

One  starts  with  the  knowledge  of  the  direction  cosines  of  OT  (x  ,  y  ,  and  z  )  and  OT.  (x  ., 

6  '  n  Jn  n  l  ni 


y  .,  and  z  .).  Thus, 
‘'ni  ni 


The 

The 


cos  2u.  =  xx.+yy.  +  zz. 
ni  n  ni  */n‘/ni  n  ni 


cos  2u  . .  =  x  .x  .  +  y  .y  .  +  z  .z 
nij  ni  nj  -mi/nj  ni  nj 

cone  TOS°  =  y  intersects  the  plane  TT^T^  along  a 
part  of  the  circle  within  the  boundary  has  length 


(8a) 

etc.  (8b) 

circle  whose  circumference  is  Zn  TH.,  . 

lb 


s  =  2tt  TH.,  - 
lb 


I 

E 

i=  1 


H.  H., 
la  lb 


(9) 
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where  H.  is  the  length  of  the  arc  outside  the  segment.  The  angle  ffn  is  then 


tn 


2ttTH.,  -  2  H.  H., 
id  .  .  la  lb 
i=l 


tn  OH 


ib 


OH. 


ib 


But 


and 


TH.,  =  OH.,  siny 
ib  ib  1 


H.  H.v  =  TH.,  (tp  .  +  (p  .,  ) 

la  ib  lb  nia  mb 


=  OH.A<p  .  +</?.,)  siny 

ib  ^nia  ^nib 


Hence, 


tn 


Zir  —  y.  (</?  .  +  (p  .. ) 

^  nia  vnib 


i  1 


sin  y 


There  are  several  situations  that  must  be  considered. 


(a)  If  y  ^  u  . 

'  ni 


(b)  If  u  .  <  y  <  v  . 

ni  r  mj 


then  ip  .  -  (p  .,  =  0 

nia  nib 

,  TF. 

..  - 1  l 

then  w  .,  =  cos  ^777 — 

mb  TH 


ib 


But 


(10) 

(H) 

(12) 

(13) 


(14) 


(15) 

(16) 


TF.  =  OT  tanu  . 
1  ni 

and 

TH.,  =  OT  tany 
ib 

Hence, 


%ib 


-1 

cos 


tan  u  . 
_ ni 

tan  y 


(17) 

(18) 


09) 


then  </?..+</>  =©  ,+0 

^nib  nja  nib  nja 

Whenever  this  is  the  case,  the  angle  </>n ^  will  add  with  the  angle  </>nja 
in  Eq.  (14)  and  it  is  unnecessary  to  determine  them  separately.  The 
value  of  the  sum  is  found  by  the  same  procedure  that  is  given  in 
Appendix  B.  The  correspondences  in  the  two  cases  are 


(c)  If  y  v  . . 
r  '  mj 


OS  -  OT 
P.  ST°  -  T.TT . 0  0 


and 


Then 


1 


0 

y 


1 

2u  . 
ni 


nib 


+  © 


nja 


2u. 


nJ 


2u  . . 

mj 


nib 


+  <P 


nja 


cos  2u  .  —  cos  2u  .  cos  2u  . 
mj  ni _ nj 

sin  2u  .  sin  2u 

ni  nj 


(20) 


(21) 
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The  value  of  angle  is  derived  in  Appendix  C  and  is  given  by 


.  2 
tan  v 


2(1—  cos  2uni)  ( 1  —  cos  2unj)  ( 1  —  cos  2unij) 

~T 


t - 2 — 

niJ  1  -  cos  2umi.  -  cos  2um.  -  cos  2u^. .  +  2  cos  2un-  cos  2un-  cos  2u n- . 


(22) 


nj 


mj 


nj 


In  summary,  the  desired  probability  distribution  is  given  by  substituting  the  value  of  *tn, 
Eq.  (14),  into  Eq.  (7)  and  is 


p<“  <  ■’  =  7T  i  0 r)  £  U  ~  ^ 

V  c/  n=l  L  i=  1 


(</?  .  +  <P  .,) 

nia  nib 


N  I 

yfe) 

f‘i  H  2  ^nia+^nib1 

II 

•  H 

II 

c 

siny  dy 


siny  dy 


(23) 


The  values  of  f  are  given  in  Fig.  3,  and  those  of  (p^  are  given  by  Eqs.  (15),  (19)  or  (21). 
Note  that  the  value  of  I  may  change  with  the  value  of  n.  The  computation  may  be  simplified  if 
n  i  of  the  segments  are  alike,  another  n^  are  alike,  and  so  on.  Then 

M 

n  =  N  (24) 

m 

m=  1 


and 


P(a  <  a)  = 


I 


£ 

s 

e 

c 


N 

2 


M  I 

4-  y  n  y '  (<p  .  +  <p  ., ) 

47r  2L/  m  ^  mia  mib 

m=l  i=l 


siny  dy 


(25) 


The  values  .  are  the  values  of  <p  .  for  the  n  -like  segments.  Similarly,  values  u  .  and  v  •  • 
mi  m  m  1111 

will  be  used  for  the  values  of  u^.  and  vn^,  respectively. 

The  limits  of  integration  need  only  be  over  the  range  of  angle  y,  where  the  integrand  is  finite. 
The  lower  limit  is  then  dependent  on  the  factor  fg/lc  and  by  reference  to  Fig.  3,  its  value  is 

y^  =  0  if  a  >  p  (26a) 

and 


y  ^  -  P  —  a  if  a  ^  0  .  (26b) 

If  y  ^a  +  p,  then  fg/*c  -  0,  or  if  y  ^  vm,  where  vm  is  the  largest  value  of  vmij  for  all  values  of 

m  or  i,  then  2  (<p  .  +  <p  )  =  2t r,  and  the  integrand  is  0.  The  upper  limit  of  y  is  then 

.  mia  miD 


yu  =  a  +  p  if  a  +  P  4  vrn 


(27a) 


or 


y  =  v  if  a  +  >  v  .  (27b) 

u  m  m 

Since  the  largest  value  angle  y  can  have  is  v  ,  the  range  of  values  of  angle  a  is  0  to  P  +  vm- 
For  a  >  P  +  vm,  P (a  <  a)  =  1. 


7 


III.  PROGRAMMING  THE  EQUATIONS  FOR  COMPUTER  CALCULATION 


The  difficulties  of  programming  the  equations  for  P(a  <  a)  occur  mainly  in  the  formulation 


I 

of  the  sum  2 
i=  1 


( (p  .  +  <p  .,) . 

rmia  mib 


A  sensible  approach  is  indicated  in  the  following  scheme. 


Having  a  value  of  y,  we  determine  whether  it  is  larger  than  each  of  the  values  of  v  ^  , , 

v  v  ^  ,,  etc.,  in  turn.  In  each  case  where  y  is  larger,  the  corresponding  value  of  </>  .,  + 

m23  m34  b  mib 

(p  .  [Eq.  (21)]  is  found  and  added.  Then  we  compare  the  value  of  y  with  each  of  the  values 

of  u  ,,  u  u  etc.  If  y  is  smaller,  there  is  no  contribution  f Eq .  (15)].  If,  for  example, 

the  value  of  y  is  larger  than  u  both  values  <p  »  and  may  contribute  to  the  sum.  If 

y  ^vm23>  the  contribution  of  (p m3a  has  already  been  taken  up  in  the  first  step.  Similarly,  if 

v^v  ,  4,  the  same  is  true  for  <p  ...  If  y  <  v  the  contribution  of  (p  ?  is  cos  \tan  u  0/tan*y) 

r  m34'  m3b  1  m23’  Y  m3a  m3'  ' 

[ Eq .  (19)].  Also,  if  y  <  vm^,  the  contribution  of  is  cos  (tan u^^/tany).  The  net  contribu¬ 

tion  then  is 


^m3a  +  ^m3b 


€  m3 


-1 

cos 


tan  u  ~ 
m3 

tan  y 


(28) 


where  e may  have  the  values  of  0,  1  or  2.  In  a  program,  one  can  start  with  the  values, 

e  .  =  2,  for  all  m  and  i.  Then,  if  it  is  found  that  y  ^  v  the  values  of  €  .  and  e  are  each 

mi  '  mij'  mi  mj 

decreased  by  1.  This  produces  the  right  contribution  for  the  sum  as  given,  for  example,  by  Eq.  (28). 

A  flow  chart  for  a  satisfactory  program  is  given  in  Appendix  D.  The  data  that  are  read  into 
the  program  include  the  pointing  directions  OT^  needed  to  define  one  segment  of  each  different 
configuration. 


IV.  SELECTION  OF  POINTING  DIRECTION  GEOMETRIES 

As  indicated  in  Sec.  I,  angle  y  has  a  maximum  value  v^  for  any  geometry,  and  the  most  de¬ 
sirable  geometry  is  the  one  that  gives  the  smallest  value  of  v  ,  permitting  the  antenna  to  be  de¬ 
signed  for  minimum  beamwidth  and  consequently  maximum  gain.  However,  no  method  is  known 
for  finding  the  optimum  geometry  for  the  general  case  of  N  antennas.  The  procedures  outlined 
in  this  section  for  various  numbers  N  may  not  give  the  optimum  geometry,  but  in  each  case,  it 
can  be  demonstrated  that  the  geometry  is  sufficiently  near  optimum  that  the  loss  in  quality  is 
small. 

For  this  study,  the  pointing  directions  will  be  represented  by  points  T^  (n  =  1,  2,  .  .  .  ,  N)  on 
the  surface  of  a  sphere  of  unit  radius,  where  they  are  actually  the  vectors  joining  the  sphere  cen¬ 
ter  to  these  points.  The  sensing  direction  will  be  represented  by  another  point  S  on  the  sphere 
surface  which  is  free  to  wander  anywhere  among  the  points  T^.  The  nearest  of  the  points  to 
point  S  is  the  one  that  defines  angle  y.  As  point  S  wanders,  it  will  see  locations  of  maximum 
y  values,  such  as  point  G^  in  Fig.  4.  Each  of  these  locations  is  spaced  by  three  equal  angles 
from  the  three  neighboring  points  of  the  set  T^.  For  if  one  travels  away  from  this  location  in 
any  direction,  angle  y  decreases.  The  N  points  may  then  be  joined  by  lines  on  the  sphere  sur¬ 
face  to  form  a  polyhedron  whose  faces  are  spherical  triangles,  and  the  number  of  vertices  is 

V  =  N  .  (29) 
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The  polyhedron  will  have  F  faces,  each  of  which  will  have  associated  with  it  a  value  of 
Vf(f  -  1,  2,  .  .  .  ,  F),  where  vf  is  the  largest  value  of  y  that  can  exist  for  the  three  points  defining 
the  triangle.  The  largest  of  the  values  vf  is  then  vm>  By  Euler's  theorem, 

F  —  E  +  V  =  2  ,  (30) 

where  E  is  the  number  of  edges.  Each  face  has  three  edges  but  each  edge  is  shared  by  two 
faces.  Hence, 

E  =  |  F  .  (31) 

By  substituting  Eq.  (29)  and  Eq.  (31)  into  Eq.  (30),  one  obtains 

F  =  2N  -  4  (32) 


and 

E  =  3N  -  6  .  (33) 

If  all  the  triangles  can  be  made  equilateral  and  congruent,  it  can  be  demonstrated  that  the 
smallest  value  of  vm  possible  will  result.  Although  there  are  only  three  polyhedra  where  this  is 
possible  (tetrahedron  N  =  4,  octahedron  N  =  6,  and  icosahedron  N  =  12),  such  a  device  is  useful 
for  any  value  of  N  because  it  permits  one  to  obtain  a  hypothetical  value  of  v  ,  say  G,  which  is 
known  to  be  smaller  than  the  best  possible  value.  If  we  have  a  geometry  that  has  a  value  of  vm 
very  close  to  the  value  of  G,  then  the  geometry  must  be  very  nearly  optimum. 

To  prove  that  the  hypothetical  geometry  is  better  than  the  optimum,  we  shall  assume  that  a 
better  geometry  than  the  hypothetical  one  exists  and  we  shall  show  that  an  impossible  situation 
follows.  For  the  assumed  geometry,  all  the  values  of  vf  would  be  less  than  the  value  G.  It  will 
be  shown  below  that  for  any  value  vf,  the  triangle  on  the  spherical  surface  with  the  largest  area 
is  equilateral.  Since  each  triangle  of  the  hypothetical  geometry  is  also  equilateral  and  vf  <  G, 
then 


where  Af  and  4ir/F  are,  respectively,  the  areas  of  a  face  of  the  assumed  geometry  and  a  face  of 
the  hypothetical  geometry.  Also,  it  is  necessary  that 


F 


L  Af  = 4,r 

f=i 


Hence, 


F 


l  ¥<4*  • 

f=l 


and  since  each  term  of  the  series  is  the  same,  then 


F(^r)  <  4ir 


which  is  impossible. 
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Fig.  5.  Geometry  showing  the  spherical  triangle  formed 
by  three  adjacent  antenna  pointing  directions  OP,  OQ  , 
and  OR. 


Consider  one  of  the  spherical  triangles,  say  PQR,  in  a  coordinate  system,  as  shown  in 
Fig.  5.  This  triangle  is  oriented  so  that  direction  OP  is  in  the  X-Z  plane  and  angles  POZ°  = 
QOZ°  =  ROZ°  =  y.  Also,  let  the  angles  separating  the  planes  containing  the  Z-axis  and  points 
P,  Q,  and  R,  be  p,  q,  and  r,  as  shown.  Hence, 

p  +  q  +  r  =  2t r  ,  (34) 

and  the  coordinates  are 

P(siny,  0,  cosy) 

Q(siny  cos  q,  siny  sinq,  cosy) 

and 


R(siny  cos  r,  —  siny  sin  r,  cosy) 

If  angles  q  and  r  have  values  where  the  area  of  triangle  PQR  is  maximum,  then  any  infinitesi¬ 
mal  change  in  q  or  r  will  produce  no  change  in  area.  That  is,  if  angle  q  increases  by  dq,  point 
Q  moves  to  Q’  and  the  change  in  area  is  0,  or 

ARQQ'  -  APQQ'  =  0  .  (3  5) 

Now  APQQ’  must  be  proportional  to  dq  to  the  first  order  and  must  be  a  function  of  q  and  y  alone. 
Hence, 


APQQ1  =  f(q,  y)  dq  (36) 

Similarly, 

ARQQ'  =  f(p,  y)  (-dp)  ,  (37) 

where  the  same  function  f  is  used  but  with  p  replacing  q.  Since  dq  -dp,  one  finds  that 

f(q,  y)  =  f(p,  y)  .  (38) 

Similarly,  a  change  in  area  due  to  an  increase  in  angle  r  is  also  0  and 

f(r,  y)  =  f(p,y)  .  (39) 

Certainly,  one  solution  of  Eqs.  (34),  (38),  and  (39)  is 

p  =  q  =  r=  y7r  =  120°  .  (40) 
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If  the  function  f  is  actually  evaluated  and  Eq.  (38)  is  solved,  one  gets  the  result 
q  =  p  +  n7r  ,  n  0,  ±1,  ±2,  ... 

However,  for  values  of  n  other  than  0,  the  area  is  smaller.  Hence,  the  triangle  with  the 
largest  area  is  equilateral. 

For  the  hypothetical  geometry,  y  =  G,  q  =  Zn/3,  and  r  =  Zir/l,  and  the  coordinates  become 
P(sin  G,  0,  cos  G) 


Q 


^3 


^  sin  G,  ^^-sinG,  cos  G 


)  ■ 


and 

R^—  ^  sinG,  -^-sinG,  cos  cj 
The  plane  containing  the  points  P,  Q,  and  O  has  equation 

x  cos  G  +  y  \T3  cos  G  —  z  sin  G  =  0  , 

and  its  normal  is 


(41) 


N  =  i  cos  G  +  j  n/~3  cos  G  —  k  sinG  .  (42) 

The  spherical  angle  QPR°  is  double  the  angle  between  N  and  the  Y-axis  and  is  given  by 


QPR  * 


j 


I  N  I 


By  substituting  Eq.  (42)  into  Eq.  (43)  and  solving  for  G,  one  obtains 

1 


_  1  4  QPR  ° 

cos  G  =  -  cot  — ^ — 


(43) 


(44) 


Since  several  lines  radiate  from  point  P  to  delineate  a  corresponding  number  of  triangles  all 
with  angle  QPR°,  there  must  be  2tt/QPR°  lines  meeting  at  each  point  P.  Referring  again  to  the 
picture  of  the  polyhedron,  we  see  that  there  are  27t/QPR°  edges  meeting  at  each  vertex  and  that 
since  each  edge  is  shared  by  2  vertices,  the  total  number  of  edges  is 


p  _  7rN 

b  "  QPR0 


(45) 


From  Eqs.  (3  3)  and  (4  5),  one  obtains 


QPR  ° 


7 r  N 

3  N  -  2 


(46) 


and  from  Eq.  (44) 

cos  G  =  — - —  cot  -r  KT  ^  7  .  (47) 

vi  6  N _2 

This  formula  provides  a  useful  method  of  indicating  how  good  the  geometry  is  that  we  choose  to 
investigate,  by  permitting  a  comparison  of  the  value  of  v^  with  that  of  G.  An  incidental  outcome 
of  this  theory  is  that  the  area  of  one  of  the  spherical  triangles  can  be  derived  in  terms  of  angle 
QPR°  by  equating  the  area  to  4tt/F  and  substituting  for  F  [Eq.  (32)]  and  then  for  N  [Eq.  (46)]. 

The  result  obtained  agrees  with  that  educed  in  spherical  trigonometry,  as  it.  should. 
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Now  we  can  proceed  to  choose  pointing  direction  geometries  and  investigate  their  character¬ 
istics  for  the  various  values  of  N.  The  values  chosen,  mainly  for  their  ease  of  solution,  are 
N  =  4,  6,  8,  12,  16,  24,  32,  48,  and  60.  In  each  case,  the  positions  of  the  points  are  chosen  in 
order  to  have  the  largest  number  of  triangles  equilateral  and  the  values  of  v^.  at  or  near  v^. 
There  is  no  evidence  that  the  choices  are  the  best,  but  the  comparisons  of  the  values  with 
G  show  them  to  be  good.  The  choices  are  based  largely  on  the  geometries  given  in  the  paper  of 
Schiitte  and  van  der  Waerden.'  In  the  case  N  =  32,  two  geometries  are  given.  One  has  slightly 
poorer  characteristics  but  its  geometry  repeats  in  the  8  octants  which  may  be  more  convenient 
for  some  applications. 


Fig.  6.  Geometry  showing  the  chosen  pointing 
directions  for  4  antennas. 


N  =  4  (Fig.  6) 

The  4  points  are  at  the  vertices  of  a  tetrahedron. 

Geometry  symmetries: 

All  triangles  are  equilateral  and  congruent. 

Coordinates: 

(0,  0,  1.00000) 

T?  (0.94281,  0,  -0.33333),  etc. 

v  =  G  =  70.53° 
m 

A  plot  of  P(a  <  a)  for  this  geometry  is  shown  in  Fig.  7. 


*  K.  Schutte  and  B.L.  van  der  Waerden,  “Auf  welcher  Kugel  haben  5,  6,  7,  8  oder  9  Punkte  mit  Mindestabstand 
Eins  Platz?"  Math.  Annalen  ^23,  96  (1951). 
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Fig.  7.  Probability  that  receiver  is  at  an  angle  less  than  a  from  the  pointing  direction 
of  the  antenna  in  operation  for  the  4-antenna  geometry  (Fig.  6). 
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/3  ( deg) 


Fig.  8.  Geometry  showing  the  chosen  pointing 
directions  for  6  antennas  . 


N  =  6  (Fig.  8) 

The  6  points  are  at  the  vertices  of  an  octahedron. 

Geometry  symmetries: 

All  triangles  are  equilateral  and  congruent. 

Coordinates: 

T  (0,  0,  1.00000) 

Tz  (1.00000,  0,  0),  etc. 

v  =  G  =  54.74°. 
m 

A  plot  of  P(n  <  a)  for  this  geometry  is  shown  in  Fig.  9. 


N  =  8  (Fig.  10) 

4  points  are  midway  between  the  X-Z  and  Y-Z  planes  and  lie  above  the 
X-Y  plane. 

4  points  are  on  the  X-Z  and  Y-Z  planes  and  lie  below  the  X-Y  plane. 
All  points  are  the  same  distance  from  the  X-Y  plane. 

Introduce  points  on  Z-axis,  on  X-Z  plane  and  on  negative 
Z-axis. 


Condition: 


T 1U 1  =  T1U2  =  T5U2  =  T5U3- 


Geometry  symmetries: 

4  triangles  (e.g.,  are  right  angle,  isosceles  and 

congruent  and  v^.  =  vm- 

8  triangles  (e.g.,  T1T2T^)  are  isosceles  and  congruent  and 

v-  =  v  . 
f  m 
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Fig.  9.  Probability  that  receiver  is  at  an  angle  less  than  a  from  the  pointing  direction 
of  the  antenna  in  operation  for  the  6-antenna  geometry  (Fig.  8). 
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Fig.  10.  Geometry  showing  the  chosen  pointing 
directions  for  8  antennas. 
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Fig.  11.  Probability  that  receiver  is  at  an  angle  less  than  a  from  the  pointing  direction 
of  the  antenna  in  operation  for  the  8-antenna  geometry  (Fig.  10). 
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Coordinates: 


T.J  (0.54120,  0.54120,  0.64359) 

T,  (0.76537,  0,  -0.64359),  etc. 

5 

U2  (0.98517,  0,  0.17157) 

v  =  49.94° 
m 

G  =  46.52°. 

A  plot  of  P(a  <  a)  for  this  geometry  is  shown  in  Fig.  11. 


N  =  12  (Fig.  12) 

The  12  points  are  at  the  vertices  of  an  icosahedron. 

Geometry  symmetries: 

All  triangles  are  equilateral  and  congruent. 

Coordinates: 

T  (0,  0,  1.00000) 

T2  (0.89443,  0,  0.44721) 

T3  (0.27639,  0.85065,  0.44721) 

T  (0.72361,  0.52573,  -0.44721),  etc. 

4 

v  -  G  =  37.38°  . 
m 

A  plot  of  P(a  <  a)  for  this  geometry  is  shown  in  Fig.  13. 

N  =  16  (Fig.  14) 

4  points  are  on  the  X-Z  and  Y-Z  planes  and  lie  above  the  X-Y  plane. 

4  points  are  midway  between  the  X-Z  and  Y-Z  planes  and  lie  below 
the  X-Y  plane. 

The  8  points  are  the  same  distance  from  the  X-Y  plane. 

4  points  are  midway  between  the  X-Z  and  Y-Z  planes  and  lie  above 
the  X-Y  plane. 

4  points  are  on  the  X-Z  and  Y-Z  planes  and  lie  below  the  X-Y  plane. 
The  8  points  are  the  same  distance  from  the  X-Y  plane. 

Introduce  points  on  Z-axis,  midway  between  X-Z  and  Y-Z 
planes  and  U^. 

Condition: 

T1U1=T1U2=T1U3=T3U2=T3U3=T4U3- 
Geometry  symmetries: 

4  triangles  (e.g.,  T^T^T^)  are  right  angle,  isosceles 
and  congruent  and  v^.  =  v  . 

8  triangles  (e.g.,  T^T^T^)  are  isosceles  and  congruent  and 

vr  =  v  . 
f  m 

16  triangles  (e.g.,  T^T^T^)  are  congruent  and  v^.  =  v^. 
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Fig.  12.  Geometry  showing  the  chosen  pointing 
directions  for  12  antennas. 


Fig.  13.  Probability  that  receiver  is  at  an  angle  less  than  a  from  the  pointing  direction 
of  the  antenna  in  operation  for  the  12-antenna  geometry  (Fig.  12). 
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Coordinates: 

T1  (0.55598,  0,  0.83119) 

T4  (0.99524,  0,  -0.09749) 

T3  (0.70374,  0.70374,  0.09749),  etc. 

U2  (0.54661,  0.54661,  0.63437) 

U3  (0.87575,  0.24798,  0.41421) 

v  =  33.78° 
m 

G  -  32.13° . 

A  plot  of  P(a  <  a)  for  this  geometry  is  shown  in  Fig.  15. 


N  24  (Fig.  16) 


3  points  are  in  each  octant  which  form  an  equilateral  triangle  centered 
in  the  octant. 

The  geometry  is  the  same  for  all  octants. 

Introduce  points  U  .  on  Z-axis,  U?  at  the  midpoint  of  the  first  octant 
and  U3. 

Condit  ions: 


T  T  T  T  -  T  T 
1  l1  2  1 2  3  13' 


T.U,  =  T,U9  =  T9U9  =  T-U,  =  T,U_  =  T7U_  -  T.U_. 
11  12  22  32  13  23  43 


Geometry  symmetries: 

12  triangles  (e.g.,  T^T^T^)  are  right  angle,  isosceles  and 
congruent  and  =  v  . 

8  triangles  (e.g.,  T,T?T3)  are  equilateral  and  congruent 

and  v-  =  v  . 
f  m 

24  triangles  (e.g.,  T1T2T4^  are  congruent  and  =  vm 


Coordinates: 


Ti  (0.26993,  0.37827,  0.88546) 

T2  (0.88546,  0.26993,  0.37827), 

U2  (0.57735,  0.57735,  0.57735) 

U3  (0.64344,  0.10754,  0.75791) 

v  =  27.69° 
m 

G  =  26.05°. 


etc. 


A  plot  of  P(n  <  a)  for  this  geometry  is  shown  in  Fig.  17. 


N  =  32,  Case  1  (Fig.  18) 

8  points  are  at  the  midpoints  of  the  octants. 

24  points  lie  on  the  coordinate  planes  at  equal  distances  from  the  axes. 
The  geometry  is  the  same  for  all  octants. 

Introduce  points  U.  on  Z-axis  and  U?  midway  between  X-Y  and  Y-Z 
plane. 
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Fig.  14.  Geometry  showing  the  chosen  pointing 
directions  for  16  antennas. 


Fig.  15.  Probability  that  receiver  is  at  an  angle  less  than  a  from  the  pointing  direction 
of  the  antenna  in  operation  for  the  16-antenna  geometry  (Fig.  14). 
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1 3- 6?- 2609 1 


Fig.  16.  Geometry  showing  the  chosen  pointing 
directions  for  24  antennas. 


Fig.  17.  Probability  that  receiver  is  at  an  angle  less  than  a  from  the  pointing  direction 
of  the  antenna  in  operation  for  the  24-antenna  geometry  (Fig.  16). 
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a  ( deg ) 


Z 


N  =  32 


l3-62-26ll| 


Fig.  18.  Geometry  for  one  set  of  pointing  directions  for  32  antennas. 


Fig.  19.  Probability  that  receiver  is  at  an  angle  less  than  a  from  the  pointing  direction 
of  the  antenna  in  operation  for  the  32-antenna  geometry  (Fig.  18). 
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Condition: 


Geometry  symmetries: 

12  triangles  (e.g.,  T^T^T^)  are  right  angle,  isosceles  and 
congruent  and  vf  =  v  . 


24  triangles  (e.g.,  T^T^T^)  are  isosceles  and  congruent  and 

vr  =  v  . 
f  m 

24  triangles  (e.g.,  T^T^T^)  are  isosceles  and  congruent  and 

vr  <  v  . 
f  m 

Coordinates: 

Ti  (0.40475,  0,  0.91443) 

T3  (0.57735,  0.57735,  0.57735) 

T4  (0.91443,  0,  0.40475),  etc. 

U2  (0.69318,  0.19747,  0.69318) 

v  -  23.88° 
m 

G  =  22.49°. 


A  plot  of  P(a  <  a)  for  this  geometry  is  shown  in  Fig.  19. 

N  =  32,  Case  2  (Fig.  20) 

12  points  are  at  the  vertices  of  an  icosahedron. 

20  points  are  at  the  midpoints  of  the  triangles  of  the  icosahedr  on. 


Geometry  symmetries: 

All  triangles  are  isosceles  and  congruent. 


1 3-62-26131 


Y 


Fig.  20.  Geometry  for  a  second  set  of  pointing  directions  for  32  antennas. 
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Fig.  21.  Geometry  showing  the  chosen  pointing  directions  for  48  antennas. 


Fig.  22.  Probability  that  receiver  is  at  an  angle  less  than  a  from  the  pointing  direction 
of  the  antenna  in  operation  for  the  48-antenna  geometry  (Fig.  21). 


24 


0  (deg) 


Coordinates: 

T 1  (0,  0,  1.00000) 

Tz  (0.89443,  0,  0.44721) 

T3  (0.49112,  0.35682,  0.79465) 

T4  (0.98225,  0,  -0.18759) 

Ts  (0.79465,  0.57735,  0.18759),  etc. 

v  -  22.69° 
m 

G  =  22.49° . 


N  =  48  (Fig.  21) 


6  points  are  in  each  octant  which  form  a  cluster  of  4  equilateral  triangles 
centered  in  the  octant. 

The  geometry  is  the  same  for  all  octants. 

Introduce  points  on  Z-axis,  at  midpoint  of  an  octant  and  U^. 


Conditions: 


T1T2  =  T2T3  =  T1T3  =  T2T4  =  T2T5  =  T4T5  =  T3T5  ;  T3T6  =  T5T6' 


T.U.  =  T,U,  =  T..U,  =  T,U,  =  T,U,  =  T,U..  =  T.U,. 

11  22  32  52  23  (3  o  3 


Geometry  symmetries: 

12  triangles  (e.g.,  T^T^T^)  are  right  angle,  isosceles  and 
congruent  and  Vj.  =  vm. 

32  triangles  (e.g.,  T1T2T3)  are  equilateral  and  congruent 
and  vf  =  vm. 

24  triangles  (e.g.,  T^T-^Tg)  are  congruent  and  vj  =  vm. 

24  triangles  (e.g.,  TjT^Ty)  are  congruent  and  Vj.  <  vm. 


Coordinates: 

(0.09252,  0.31372,  0.94500) 

T4  (0.94500,  0.09252,  0.31372) 

T,  (0.61706,  0.28702,  0.73270) 

T.  (0.73270,  0.61706,  0.28702),  etc. 

3 

U2  (0.57735,  0.57735,  0.57735) 

U3  (0.60477,  -0.03842,  0.79547) 

v  =  19.09° 
m 

G  =  18.30°  . 

A  plot  of  P (a  <  a)  for  this  geometry  is  shown  in  Fig.  22. 
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Fig.  23.  Geometry  showing  the  chosen  pointing  directions  for  60  antennas. 


Fig.  24.  Probability  that  receiver  is  at  an  angle  less  than  a  from  the  pointing  direction 
of  the  antenna  in  operation  for  the  60-antenna  geometry  (Fig.  23). 
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/3  (deg) 


N  -  60  (Fig.  23) 


36  points  lie  on  the  coordinate  planes. 

24  of  the  36  points  are  the  same  distance  from  the  axes. 

The  other  12  are  midway  between  2  axes. 

24  points  are  midway  between  two  coordinate  planes  at  equal  distances 
from  the  centers  of  the  octants. 

The  geometry  is  the  same  for  all  octants. 

Introduce  points  U  .  on  X-axis,  U?  midway  between  X-Y  and  Y-Z  planes 
and  U3. 

Condition: 

T3U1  =  T3U3=T2U3  =  T11U3  =  T2U2  =  T11U2- 
Geometry  symmetries: 

12  triangles  (e.g.,  T^T^T^)  are  right  angle,  isosceles 
and  congruent  and  Vf  =  v^. 

24  triangles  (e.g.,  T^T^T,,)  are  isosceles  and  congruent 
andvf<vm. 

48  triangles  (e.g.,  T^T^T^)  are  congruent  and  vf  =  vm- 

24  triangles  (e.g.,  T^T^T^)  are  isosceles  and  congruent 

and  vr  =  v 

f  rn 

8  triangles  (e.g.,  TioTllTiz^  are  e(4uilateral  and  congruent 
and  Vf  <  v  . 

Coordinates: 

T3  (0.95660,  0,  0.  29140) 

Tz  (0.70711,  0,  0.70711) 

T  .  .  (0.77565,  0.44630,  0.44630),  etc. 

1 1 

U2  (0.67642,  0.29140,  0.67642) 

U3  (0.84544,  0.16667,  0.50740) 

v  =  16.94° 
m 

G  =  16.35° . 

A  plot  of  P(a  <  a)  for  this  geometry  is  shown  in  Fig.  24. 

V.  INCLUSION  OF  ANTENNA  PATTERN  CHARACTERISTICS 

Having  determined  the  probability  distribution  of  angle  a  for  a  specific  geometry,  one  will 
want  to  design  an  antenna  whose  pattern  is  most  suitable  for  the  case  of  interest.  As  an  example, 
suppose  we  want  to  build  a  satellite  with  32  antennas  arranged  as  shown  in  Fig.  18,  which  is  to 
have  an  RF  power  source  strong  enough  to  allow  messages  to  be  received  100  percent  of  the  time 
at  stations  located  anywhere  where  0  is  8°  or  less.  Referring  to  Fig.  19,  we  see  that  angle  a 
will  have  values  up  to  a  -  31.9°.  The  antenna  is  designed  for  maximum  gain  at  this  angle,  thus 
minimizing  the  power  requirements  of  the  RF  source.  If  reception  for  80  percent  of  the  time 
were  allowed,  angle  a  would  be  less  than  a  =  21.6°  for  80  percent  of  the  time,  and  the  antennas 
could  be  designed  for  maximum  gain  at  21.6°  which  would  be  higher  than  the  gain  of  the  first 
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antenna.  Similarly,  smaller  values  of  /3  would  permit  smaller  values  of  a  and  correspondingly 
higher  antenna  gains  for  the  same  probability. 

To  illustrate  better  the  application  of  the  antenna  pattern,  suppose  the  pattern  is  given  by 

2<J  i ( gQ  sin  a) 


where  J  ^  is  the  Bessel  function  of  the  first  kind,  g  is  the  signal  amplitude  at  angle  a  off  the 
axis,  and  gQ  is  the  signal  amplitude  on  the  axis  (a  -  0) .  This  pattern  is  based  on  the  situation 
where  a  plane  wavefront  emanates  from  a  circular  aperture  and  it  is  fairly  representative  of 
RF  horns  in  general.  The  antenna  designer  has  a  choice  of  the  value  gQ  at  his  disposal.  It  is 
approximately 


Ztt  r 
o  “  X 


(49) 


where  r  is  the  aperture  radius. 

The  half-power  beamwidth  HPBW  is  found  by  Eq.  (48),  where 

a  =  \  HPBW  (50) 


and 


(51) 


ANGLE  a  (deg) 


Fig.  25.  Pattern  of  a  circular  pipe  of  1-inch  aperture  butted  against  a  Luneburg  lens 
of  3- inch  diameter  taken  at  a  frequency  of  7750  Mcps. 
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Hence, 


or 


g  sin  4-  HPBW  =  1.616 
6o  2 

HPBW  =  2 


(52) 


(53) 


are. 


and 


Now  at  a  -  a,  gain  g  is  maximum  when  gQ  has  a  value  where  J  ^  is  maximum, 
from  Eq.  (48), 

g  sin  a  =  1.841 
6o 

g  sin  a  =  1. 164 


The  values 

(54) 

(55) 


An  immediate  observation  is  that 


JL  =  0.632  -  —  3.98  db  ,  (56) 

go 

and  the  maximum  range  of  gain  values  is  3.98  db  for  the  optimum  pattern.  In  our  example, 
where  a  =  31.9°,  gQ  -  3.486  -*•  10.85db,  and  g  =  2.204  -*  6.87 db. 

A  similar  analysis  was  applied  to  the  other  geometries  and  the  results  are  summarized  in 
Appendix  A.  Two  values  of  /?,  0°  and  8°,  were  considered.  The  trends  in  a  change  of  N  or  (3 
are  easily  noted  and  for  values  not  included,  the  characteristics  can  be  estimated. 

It  is  useful  to  give  an  example  of  the  application  of  an  actual  pattern  to  a  probability  distri¬ 
bution.  The  pattern  shown  in  Fig.  25  (very  nearly  the  same  in  the  E-  or  H-plane)  is  of  a  circular 
pipe  of  1-  inch  aperture  butted  against  a  Luneburg  lens  of  3 -inch  diameter  taken  at  a  frequency 
of  7750  Mcps.  For  this  pattern,  HPBW  =  34°  and  gQ  =  5.69  -**  15.1db,  and  hence  gQ  sin|  HPBW 
1.66,  which  agrees  well  with  the  value  given  by  Eq.  (52).  This  pattern,  in  conjunction  with  the 
geometry  of  Fig.  18,  produces  a  probability  distribution  as  shown  in  Fig.  26.  These  curves  were 


Fig.  26.  Probability  that  the  antenna  gain  is  less  than  a  value  g*  in  the  direction 
of  the  receiver  for  the  geometry  of  Fig.  18  using  the  pattern  in  Fig.  25. 
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|  3  -  6?- ?  6  ?Q(o  -  eT] 
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Elements 
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SOURCE 
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4 

0.6 

0.  1 

ANTENNAS 

4 

0.8 

RF  SOURCE 


Switching  Estimated 
Elements  Loss  (db) 

0.  1 

2  0.3 

0.  1 

6  0.45 

0.  1 
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(a)  4-antenna  system. 


(b)  6-antenna  system. 


r  f  source 


Switching 

Elements 


10 


Estimated 
Loss  (db) 

0.  1 

0.3 

0.  1 

0.6 

0.  1 


(c)  8-antenna  system. 


1.2 


RF  SOURCE 


Switching  Estimated 
Elements  Loss  (db) 

0.  1 

3  0.45 

0.  1 

12  0.6 

0.  1 


15 


(d)  1 2-antenna  system  . 


1.35 


RF  SOURCE 


Switching 

Elements 


4 


16 


20 


Estimated 
Loss  (db) 

0.  1 
0.6 

0.  1 

0.6 

0.1 

1.5 


(e)  16-antenna  system. 


Fig.  27(a-i).  Switching  matrix  and  estimated  insertion  loss  for  various  values  of  N. 
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simply  obtained  by  converting  the  angles  a  of  Fig.  19  to  gains  by  Fig.  25.  The  choice  of  pattern 
for  this  geometry  is  not  optimum,  and  where  the  gains  may  drop  to  say  8.5  db  for  (3  =  0°  or 
2.1db  for  (3  =  8°,  the  gains  would  drop  at  the  most  to  9.2  and  6.9  db,  respectively,  if  the  optimum 
patterns  were  used. 

The  ideal  antenna  for  this  application  is  one  where  the  gain  is  constant  for  angle  a  less  than 
a  and  0  otherwise  (within  the  restriction  of  retaining  circular  symmetry).  Then 

^ideal  •  JL  *  n  <  a  (57a) 


or 


^ideal 


n  >  a 


(57b) 


On  comparison  with  the  gain  given  by  Eq.  (55), 


a 

g.  ,  ,  COS  T 

foideal  _  2 

g  “  0.582 


(58) 


Approximately,  cos  a/2  =  1  and  the  ratio  is  4.7  db.  If  such  a  pattern  could  be  created,  the  RF 
power  necessary  for  communication  is  correspondingly  reduced.  However,  the  opportunities 
for  shaping  a  more  desirable  pattern  are  impeded  by  the  difficulties  of  maintaining  a  simple 
antenna  design. 


VI.  CONSIDERATION  OF  THE  RF  CIRCUITS 

To  complete  this  study,  some  consideration  should  be  given  to  the  switching  matrix  supply¬ 
ing  the  RF  to  the  appropriate  antenna.  The  matrix  will  be  assumed  to  contain  switches  that  are 
either  SP2T,  SP3T  or  SP4T.  The  expected  switching  matrices  for  the  various  values  of  N  are 
shown  in  Fig.  27(a-i).  If  we  assume  that  there  is  an  insertion  loss  of  0.6  db  in  a  SP4T  switch, 

0.45  db  in  a  SP3T  switch,  0.3  db  in  a  SP2T  switch  and  0.1  db  in  each  length  of  transmission  line 
connecting  the  switches  to  each  other,  to  the  power  source  or  to  the  antennas,  we  can  estimate 
the  total  loss  of  RF  through  the  matrix.  The  results  are  given  in  Appendix  A.  (The  values 
of  switch  insertion  losses  are  based  on  experimental  evidence  obtained  for  switches  using  PIN 
diodes  in  waveguide  packages,  where  f  =  7750  Mcps.  Ferrites  should  show  better  characteristics 
but  they  are  less  desirable  for  satellite  applications  due  to  weight  and  driving  power  requirements.) 

Each  SP2T,  SP3T  or  SP4T  switch  requires  2,  3  or  4  switching  elements.  For  example,  for 
the  matrix  of  32  antennas,  42  switching  elements  are  needed.  One  can  appreciate  the  complexity 
of  the  logic  circuits  and  the  problems  of  generating  power  to  operate  switching  elements,  if  the 
number  of  antennas  is  large. 

VII.  CONCLUSIONS 

The  results  of  this  study  are  well  summarized  in  Appendix  A.  The  choice  of  the  number  of 
antennas  will  also  depend  on  factors  other  than  those  mentioned  in  this  report.  However,  it  ap¬ 
pears  that  for  the  small  reduction  in  RF  power  requirements,  more  than  12  or  16  antennas  are 
not  warranted,  due  to  the  complexity  of  the  switching  circuits.  A  fair  appraisal  should  also  in¬ 
clude  a  comparison  of  this  type  of  satellite  with  one  that  does  not  contain  switches  but  is  de¬ 
signed  to  radiate  as  close  to  isotropic  as  possible.  In  any  case,  this  system  is  put  forth  as  a 
preferable  means,  for  some  types  of  communication,  of  linking  two  distant  locations  on  the  earth. 


31 


3-6?-26?0(f -  g)  | 


RF  SOURCE 


Switching 

Elements 


Estimated 
Loss  (db) 

0.  1 

0.3 

0.  1 

0.45 

0.  1 

0.6 
0.  1 

1.75 


(0  24  -antenna  system. 


Switching 

Elements 

RF  SOURCE 


Estimated 
Loss  (db) 

0.  1 

0.3 

0.  1 

0.6 

0.  1 
0.6 
0.  1 
1.9 


(g)  32-antenna  system. 


Fig.  27.  Continued. 
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Fig.  27.  Continued . 
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APPENDIX  A 


SUMMARY  OF  EXPECTED  CHARACTERISTICS  FOR  CASES  INVESTIGATED 


Number  of  Antennas  N 

4 

6 

8 

12 

16 

24 

32 

Case  2 

32 

Case  1 

48 

60 

Number  of  switching  elements 

4 

8 

10 

15 

20 

32 

42 

42 

63 

79 

Estimated  circuit  loss  (db) 

0.8 

1.0 

1.2 

1.3 

1.5 

1.8 

1.9 

1.9 

2.0 

2.2 

P  =  0°  (angle  between  receiver  and  sensing 
direction) 

Largest  possible  angle  of  receiver 
from  axis  of  operating  antenna 
(deg) 

70.5 

54.7 

49.9 

37.4 

33.8 

27.7 

22.7 

23.9 

19.  1 

16.9 

This  angle  if  angles  between 
adjacent  pointing  directions 
could  be  made  equal  (ideal 
case)  (deg) 

70.5 

54.7 

46.5 

37.4 

32.  1 

26.  1 

22.5 

22.5 

18.3 

16.3 

Maximum  antenna  gain  at  this 
angle  (db) 

1.83 

3.08 

3.65 

5.66 

6.42 

7. 98 

9.60 

9.18 

11.03 

12.  03 

Gain  at  this  angle  after 
deducting  circuit  loss  (db) 

1.0 

2.1 

2.4 

4.4 

4.9 

6.2 

7.7 

7.3 

9.0 

9.8 

Antenna  gain  at  peak  for  same 
antenna  pattern  (db) 

5.81 

7.06 

7.63 

9.64 

10.40 

1 1.96 

13.58 

13.16 

15.01 

16.01 

Gain  at  peak  after  deducting 
circuit  loss  (db) 

5.0 

6.  1 

6.4 

8.3 

8.9 

10.2 

11.7 

11.3 

13.0 

13.8 

Half-power  beamwidth  for 
same  pattern  (deg) 

111.7 

91.6 

84.4 

64.4 

58.4 

48.2 

39.6 

41.6 

33.4 

29.6 

P  =  8° 

Largest  possible  angle  of  receiver 
from  axis  of  operating  antenna 
(deg) 

78.5 

62.7 

57.9 

45.4 

41.8 

35.7 

30.7 

31.9 

27.1 

24.9 

Maximum  antenna  gain  at  this 
angle  (db) 

1.50 

2.  35 

2.  76 

4.27 

4.85 

6.00 

7.  16 

6.87 

8.  15 

8.  82 

Gain  at  this  angle  after 
deducting  circuit  loss  (db) 

0.7 

1.3 

1.6 

3.0 

3.3 

4.2 

5.3 

5.0 

6.  1 

6.6 

Antenna  gain  at  peak  for  same 
antenna  pattern  (db) 

5.48 

6.  33 

6.74 

8.  25 

8.83 

9.98 

11.14 

10.  85 

12.  13 

12.  80 

Gain  at  peak  after  deducting 
circuit  loss  (db) 

4.7 

5.3 

5.5 

6.9 

7.3 

8.2 

9.2 

8.9 

10.1 

10.6 

Half-power  beamwidth  for 
same  pattern  (deg) 

118.7 

102.6 

96.  1 

77.3 

71.6 

61.6 

53.2 

55.2 

47.  1 

43.5 
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APPENDIX  B 

DERIVATION  OF  THE  FORMULA  FOR  f  /f  WHERE  |a  -  j3|  <y<a  +  /3 


If  a  —  /3  <  y  <  a  +  p  and  a  >  p ,  or/3  —  a<y<a+/3  and  a  p ,  then  there  will  be  two  directions 
OP1  and  OP2  on  the  cone  ROS°  -  p  where  POT°  -  a  (Fig.  B-l).  Let  points  Pjt  S,  and  T 

lie  on  a  plane  perpendicular  to  OS.  The  ratio  I  /lc  becomes  simply  the  ratio  of  the  angle  sub¬ 
tended  by  the  arc  P^P2  at  S  to  Zn ,  or 

£  2P .  ST° 

si  *  i~>  j  \ 


Now 

pTs2  +  st2  -  pTt2 

P.ST°  =  cos"  — - == - =r^ -  .  (B-2) 

ZP^  •  ST 

Also 

P^T2  =  OP^2  +  OT 2  -  20P^  •  OT  cos  a  .  (B-3) 


Hence, 

,  PTs2  +  ST2  -  OP.2  -  OT2  +  20PT  •  OT  cos  a 

n  __  -1  1  1 


P  ST°  =  cos'  — - - - - — - - -  (B-4) 

2P.S  •  ST 
1 

But 

P^S  =  OS  tan/3  ,  (B-5) 

ST  =  OS  tan  y  ,  (B-6) 

OPj  =  OS  sec  p  (B-7) 

and 

OT  -  OS  secy  .  (B-8) 
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Hence, 


P  .STC 


-  cos 


Zn 
tan  (3 


2  2  2 

tan  y  ~~  sec  P  ~~  sec  7  +  2  sec  (3  sec  y  cos  a 
2  tan  p  tan  y 


and  this  simplifies 


Therefore, 


P  ^  ST 0 


-1  .  cos  a  —  cos  P  cos  y, 

cos  ( - r — -= — t— ~ - r-) 

sin  j8  sin  y 


f 

s 

f 

c 


1  -1  cos  a  —  cos  B  cos  y 

—  COS  - i - zz - r— - - 

7r  sin  p  sin  y 


(B-9) 


( B  -10) 
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APPENDIX  C 

DERIVATION  OF  THE  VALUE  OF  THE  ANGLE  v  .. 


Referring  to  Fig.  4,  one  begins  with  the  relation 


cos  © 


nib 


TF. 

1 

TO. 

1J 


Hut 


TFi  =  OT  tanUni 


and 


Hence , 


Similarly, 


TG..  =  OT  tan v  .. 
ij  nij 

tan  u  . 
„  ni 

cos  © 

nib  tanv  . 

niJ 


cos  © 


nja 


tan  u  . 

_ nj_ 

tanv  .. 
mj 


Take  the  identity. 


sin©  sin©  -cos©  ..  cos  ©M.  -  cos  (©  .,  +  ©  .  ) 

nib  nja  nib  nja  nib  nja 


square,  substitute 


2  2 
sin  ©  1  —  cos  ©  . 

mh  m 


nib 


nib 


and 


2 

sin  © 


nja 


i  —  cos  ©  . 

nja 


and  simplify.  The  result  is 


1  —  cos^  ©  —  cos^  ©  -  —  2  cos  ©  .,  cos  © 

m  h  nj  a  n  i  h 


COS  (©  .,  +  ©  ) 

nja  nib  nja 

+  cos^  (©..+©.) 

'  nib  nja 

When  we  substitute  Eqs.  (C-4)  and  (C-5)  and  solve  for  tan^  vn[j*  ^ie  result  is 

2  2 

0  tan  u  .  +  tan  u  .  —  2  tanu  .  tanu  cos  (0  -it  ©  •  ) 

2  ni  nj  ni  nj  nib  nja 

tan  v  ..  =  - - - 5 - - - 


nij 


1  —  cos^  (©  ..  +  ©  ) 


But  from  Eqs.  (20)  and  (21) 


c°s  (enib 


©  .  )  = 
nja' 


cos  2  u  . .  —  cos  2  u  .  cos  2  u 
_ mj _ ru _ nj 

sin  2  u  .  sin  2  u 

ni  nj 


(C-1) 


(1-2) 


(C-  3) 


(C-4) 


(C-5) 


(C-6) 


(C-7) 


(C-8) 


(C-9) 


(C-10) 


(C-ll) 
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Also 


tan  u 


1  —  cos  2  u  . 

_ ni 

sin  2  u  . 
ni 


(C-12) 


and 


tan  u 


1  —  cos  2  u 

JR 


nj 


sin  2  u 


nJ 


(C-13) 


Equations  (C-ll)  to  (C-13)  are  substituted  into  Eq.  (C-10),  the  numerator  and  denominator  are 
2  2 

multiplied  by  sin  2  u^.  sin  2  u^. ,  and  the  result  is 

tan2  v  •  •  —  (1  cos  2u  .)2  sin^  2u  t  (1  -  cos  2u  )2  sin2  2u  . 
nij  ni  nj  nj  ni 


—  2  ( 1  —  cos  2u  . )  ( 1  —  cos  2u  )  (cos  2u  . .  -  cos  2u  .  cos  2u  ) 
_ m _ nj _ try _ ni _ nj_^ 

sin2  2u  .  sin2  2u  —  (cos  2u  . .  —  cos  2u  .  cos  2u  )2 

m  nj  nij  ni  nj 


(C-14) 


By  substituting 


.  2  _  ,  2  , 
sin  2u  .  1  —  cos  2u  . 

ni  ni 


(C-15) 


and 


.2  -  4  2  . 
sin  2u  1  —  cos  2u 
nj  nj 


(C-l  6) 


Eq.  (C-14)  simplifies  to 


tan  v  . 
mj 


2(1-  cos  2u  .)  1  -  cos  2u  )  1  —  cos  2u  .  ) 
_ ni  _ nj _ niT _ 

2  2  2 

1  -  cos  2u  .  -  cos  2u  -  cos  2u  . .  4  2  cos  2u  .  cos  2u  cos  2u  . 

ni  nj  nij  m  nj  nij 


.  (C-l  7) 
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APPENDIX  D 

COMPUTER  PROGRAM  TO  CALCULATE  PROBABILITY 
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Reading  input  data 


(o)2292-29-£ 


$  E 


a  u 
E  a 
•"  a> 
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Setting  of  value  of  j 


(q) 2292-29-f 
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Integration  loop 


(P)2292-29-£ 


/ 
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